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'"T. ABSTRACT 

, Aims. We use the VVV survey observations in bulge regions close to the Galactic plane to trace 

(~| the bar inclination at the Galactic latitude b ~ ±1 and to investigate a distinct structure in the 

inner regions of the bar that was previously detected at positive latitude (b = +1). 
O ' Methods. We use the (J-K^) colors of the red clump stars to obtain reddening values on 6 x 6 ar- 

' cmin scale, minimizing the problems arising from differential extinction. Dereddened magnitudes 

are then used to build the luminosity function of the bulge in regions of ~0.4 sq deg to obtain the 
mean red clump magnitudes. These are used as distance indicators to trace the bar structure. 
Results. The luminosity function clearly shows the red clump mean magnitude variation with 
1/^ , longitude, as expected from a large scale bar oriented towards us at positive Galactic longitude, 

^ with a dereddened magnitude varying from K^g = 13.4 at / = -10° to K^g = 12.4 at / = +10°. 

' We detect a change in the orientation of the bar in the central regions with |/| < 4° at b = ±1°, 

in agreement with results obtained at positive latitudes by other authors. Our results are based on 
a different dataset and at different latitude, which shows that this change in the bar orientation is 
real. This suggests that there is an inner structure distinct to the large-scale Galactic bar, with a 
different orientation angle. This inner structure could be a secondary, inner bar, with a semi-major 
, axis of ~ 500 pc that is symmetric with respect to the Galactic plane. 

H ■ 

' Key words. Galaxy: structure - Galaxy: bulge 

1. Introduction 



> 



The presence of bars in spiral galaxies is common (lEskridge et al.ll2000l) and the Milky Way is no 
exception. Different studies based on various techniques have suggested that there is a large-scale 
bar \ yith ~2.5 kpc in radius oriented 15-30 degrees with respect to the Galactic center line of sight 



DweketalJll995 



Cabrera-Lavers et al 



Stanek et al 



2007 



1997 



Bissantz & Gerhar dbooj Babusiaux & GilmorjboOS 
T — 

Rattenburv et alJI2007ft with the near end towards positive longitudes. 
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Fig. 1. K-band extinction map in Galactic coordinates showing the fields analyzed in this study. 
The 28 VVV tiles cover ranges of longitude -10° < / < 10° and of latitude -1.4° <b < -0.4° and 
0.7° <b< 1.7°. 



However, observations of external galaxies have led to th e conclusion that the existence of 



Laine et al. 



2002; 



Erwin 



201 11) and thus have provided important 



double bars is not unusual (e.g 
constraints on their scales. Secondary bars are found to be gene rally srnall, w ith a typical size of 
~ 500 pc or a relative length of about 12% of that of the main bar ( Erwinl201 1 ). Some observations 
in the very inner regions of the Milky Way hav e indeed suggested that there is an o ther structure 



Alard2001 



Nishivama et al 



2005) 



AlardlOOOlh 



inside the Galactic bar, namely a nuclear bar (e.g. 
used 2MASS star counts to show evidence of an excess in the inner regions of the Milky Way, which 



Nishivama et al 



(2005 



is interpreted as the presence of a small, lopsided, nuclear bar. In particular, 
hereafter N05) used the magnitude of red clump (RC) stars at a latitude of b = +1° to measure the 
properties of the inner regions (-10° < I < 10°) detecting a shallower inclination angle for the 
Galactic bar at fields with longitudes |/| < 4°. This was interpreted by N05 as the effect of a distinct 
inner bar structure. The existen ce of such a nuclear bar has strong effec ts on the gas distribution 
in the inner parts of the Galaxy. iRodriguez-Fernandez & CombesI (120081) concluded that indeed a 
nuclear bar is likely to be responsible for the observed so-called central molecular zone. 

In this Letter, we make use of the ESO Vista Variables in the Via Lactea survey (VVV) to trace 
the magnitudes of RC giants in the longitude range -10° < / < 10° at = ±1° to constrain the mor- 
phology of the Galactic bar. The technique, based on using RC giants a s accurate distance indicators 

(Il997h and adopted later by several 



Stanek et al 



across different lines of sight, was first applied by 
authors to trace the bar morphology at different Gal a ctic latitudes (e.g. 



Rattenburv et al. 


2007 


Cabrera-Lavers et al. 


2007) 



Minniti et al 



Babusiaux & Gilmor el2005l 



i 



(1201 11) applied this technique to 



the VVV survey data to detect the edge of the Galactic stellar disk. Here, we use it to test whether 
the change in the slope in the orientation of the Galactic bar observed by N05 at b = -1-1° is also 
evident in our data at b = -t-l°, as well as at b = -1°. We corroborate the measurements of N05, 
which suggests that there is a distinct structure in the inner regions of the Galactic bar that appears 
to be symmetric with respect to the Galactic plane. 



2. The data 

For this study, we use the near-lR JK photometry fr om the VVV public surve}U. A detailed de- 

(I2OIOI) . In part icular, the pho- 



Minniti et al 



scription of the survey and the data can be found in 
tometric calibration and dereddened magnitudes were calculated as described in iGonzalez et al 



http ://v V vsurvey. org 
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(1201 ih . hence are only briefly described here. Multi-band catalogs were produced by a cross- 
match of sources between single-band catalogs produ ced at the Cambridge Astronomical Survey 
Unit (CASU) using the STILTS code jTavloJbooel) . For this work, we analyzed a total of 28 
tiles, containing about 30 million measured sources, which cover the survey region between 
-1.4° < b < -0.4° and +0.7° < b < +1.1° at longitudes across the range -10° < Z < 10° 
(FigID- The photometric calibration was obtained by cross-matching sources of VVV with those 
of 2MASS flagged with high quality photometry. This produced the final VVV J, H, and Ks pho- 
tometric catalogs fully consistent with the 2MASS photometric system. The calibrated VVV tile 
catalogues were then dereddened and used to build the luminosity functions for each of the tiles. 



3. Analysis and results 



The extinction correction was done following the pr escriptions described inlGonzalez et al 



However, for the present analysis, we adopted the iNishivama et al 



Cardelli et al 
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(1201 Ih . 



(2009h extinction law, where 



(Il989l) . as the latter does not 



At - 0.528 • E( J- Ks) instead of the standard values of|( 
seem to be consistent with observations in the high-reddening regions analyzed here. Therefore, 
the dereddened Ks„ magnitudes were obtained as 



Ks„ ^Ks + 0.528[(7 - Ks)q -(J- K^)], 



(1) 



where w e use the mean intrinsic RC color for the bulge of (J - Ks)o - 0.68 as measured in Baade's 
window dGonzalez et al.ll201 lb . We note that by following this procedure we obtain reddening 
values in subfields of 6 x 6 arcmin instead of assigning individual extinctio n values to each star 
as in N05. Owing to the broadness of the bulge metallicity distribution (e.g. IZoccali et alj|2008t 
Bensbv et al a single intrinsic RC color cannot be adopted for all stars but only to the mean 



of the observed RC color distribution in small subfields where the effects of dififerential extinction 
are minimized. Therefore, our procedure allows us to obtain dereddened Ks magnitudes while 
keeping the imprints of the original bulge stellar population. Figure [1] shows the final extinction 
map and the field coverage for this study. 

We applied a color cut [J - Ks)o > 0.4 to select the region of the color-magnitude diagram 
(CMD) dominated by bulge giants and avoid contamination from disk stars. Using the dereddened 
Ksg magnitudes of bul ge giants, we then bu ilt the luminosity function for each tile and followed the 
method introduced by 



Stanek et al 



(Il997h to measure the RC mean magnitudes in order to trace 
the mean distance of the bulge population. The base of the luminosity function was fitted with a 
second order polynomial and a Gaussian fit was then appUed to measure the RC mean magnitude. 
Therefore the final fitting function is 

Nrc r (Ki^-Ksof 



NiKs,) = + bKs, + cK' + 



: exp 



(2) 



o-Rc V27r 

The results of the Gaussian fits for each RC are listed in Table [T] We also included an additional 
Gaussian fit for a peak in the luminosity function, distinct from the RC, which is centered at Ks^ ~ 
13.7. N05 mention the detection of this additional peak at a dereddened magnitude Kq ~13.5 but 
they only perform ed Gaussian fits up to a magnitude where the RC is unaffected by this peak. 



Natafetal 



(1201 11) presented its detection, b ased on optical OGLE data, and identified it as an 



anomalous red giant branch bump (RGBb). In 



Gonzalez et al 



(1201 ih . the peak was observed using 
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Fig. 2. Ksf) luminosity functions for representative fields at = - 1 ° (left panels) and h - (right 
panels). The red solid line shows the best fit to each distribution while the dashed line shows the 
individual Gaussian fits to the RC and to the additional peak when detected. 



the VVV photometry up to = -4°. For larger d istances from the plane, the d etection of this peak 
becomes difficult because of the split in the RC (iMc William & ZoccaUlbOlOh . However, we point 
out that, while the magnitude of the RC becomes brighter toward positive longitudes according to 
the bar orientation, the Ks^ magnitude of this additional peak changes in the opposite direction. 
This is not expected if this feature corresponds to the RGBb. A possibility could be that it instead 
corresponds to RC stars from a more distant population. Assuming an intrinsic RC magnitude 
of M/f - -1.55, the observed fainter RC peak would be at ~ 11.2 kpc from the Sun. Further 
discussion of this point, such as the number counts of the RC with respect to those of the secondary 
peak cannot be addressed in this work because of the possible incompleteness at the magnitudes 
where this peak is observed, in particular for regions affected by high reddening and crowding. 
However, we assume here that the RC mean magnitudes remain reliable as long as the additional 
peak is included in the fitting procedure, independently of its origin. Figure |2] shows the results 
for the fits to the ^T^ o distribution of RC stars in representative fields of the regions studied in this 
work. 
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Figure |3] shows the measured location of the Galactic bar obtained from the dereddened Ks„ 



magnitude of the RC alb - + 1 ° , compared with those of N05 for b_ 
magnitude of the RC of = -1.55 following 



Gonzalez et al 



- 1 ° and adopting an intrinsic 



(1201 11) . Along every line of sight, 
the dispersion in distanc e was obtained from the mea sured cr corrected by an intrinsic dispersion 
in the RC o f 0.17 mag (iB a busiaux & Gilmord l2005h and a photometric error of 0.05 mag. As 
discussed in 



Stanek et al 



119941) . the observed orientation angle based on our method differs from 
the real one by an amount that depends on the thickness of the bar, measuring larger angles for a 
thicker bar. For this reason, we included in Figure [3]examp les of the observed orientations of a bar 



Stanek et al. 



JSe 



with axis ratios of x/y = 0.15 (iCabrera-Lavers et alJl2007h for the true angles of 15°, 30°, and 45 
following corrections shown in 



(119941) . Results are in good agreement with previous 
studies, which are consistent with the bar being located at 7.6 kpc, oriented at 30°, and with its 
near end towards positive longitudes. However, the position angle of the bar in the inner ~1 kpc is 
more than 20° larger than measured in the outer regions, which implies that there is an inner central 
structure. 



4. Discussion 



Recent studies seem to point out that secondary bars in barred spirals are inde ed a common phe 



Laine et al 



2002). 



nomenon, with at least 29% of baiTed galaxies having a secondary bar (e.g 
Therefore, it would not be surprising that our Milky Way hosts such an inner structure. In this 
work, we have presented further evidence of the detection of an inner structure in the Galactic bar, 
as previously suggested by N05, identified by a change in the orientation angle of the bar measured 
at both positive (b - +1) and negative {b = -1°) latitudes, since our results have been obtained 
independently and based on a different dataset from that of N05, we confirm that the detection is 
not an artifact owing to possible issues with either the photometry or calibrations. 

However, we note that both studies are based on the use of the observed magnitude of the RC to 
trace distances across different regions of the bar. Th erefore, one point to con sider is the dependence 
of metallicity for the intrinsic magnitude of the RC ( Girardi & Salarij2001 ). This dependence can 
produce changes of ~0.1 mag for 0.2 dex changes in [Fe/H]. H owever, these gradie nts in [Fe/H] 



20081) and are even 



have only been observed along the minor axis in the outer bulg e dZoccali et al. 
more likely to be absent from the inner regions around b 1 ° ( Rich et alj2007l) . Even if gradients 



are indeed present also along the major axis for the inner regions considered here, to produce a 
similar effect to what we see in Fig. [3] the [Fe/H] gradient could not be radial, but should go from 
one extreme of the bar to the other with a slope that changes in the inner longitudes \l\ < 4°. These 
particular properties for the stellar populations of the Milky Way bar are very unlikely. 

A second point to consider is the adopted extinction law for the determination of dereddened 
Ksf, magnitudes. For a consistent comparison with N05 measurements, we have used the same 



extinction law based on 



Nishivama et al 



(120091) . However, we note that this extinction law, which 



Cardelli et al 



1989 



Savage & Mathis 



1979 



differs from the more co mmonly adopted laws (e.g. 

Rieke & Lebofsk vl l 19851 was derived for the highly reddened inner bulge fields {\b\ < 2°), while 
this does not seem to be the case for the whole bulge. In the study of RR Lyrae in the bulge 
from 



Kunder et al. 



(I2008h . variations with respect to the standard reddening law were observed 
in a few directions. However, they conclude that, on average, the extinction law of Ry ~ 3.1 is 
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Table 1. Fitting parameters for the RC and the additional peak 
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Notes. Listed are the mean dereddened (A'^^), sigma (cr) and peak of the Gaussian used to fit the RC in each 
one of the tiles. K^^^^, peak2 and 0-2 are the parameters used fo fit the secondary peak. 



consistent with the observations of the bulge at larger distances from the Galactic plane. Given this 
uncertain ty in the correct extinct ion law with line of sight, a variation of it from standard values 
to that of iNishivama et al.l (l2009h could produce a change of up to 0.1 magnitudes in /T^^ for the 
regions with the highest reddening (A^ ~ 2.5), as seen in Fig. [T] However, these variations are 
Ukely to be randomly distributed along different lines of sight. That the same change in the bar 
orientation is now observed at both positive and negative latitudes, where reddening patterns differ 
greatly (Fig.[TJ, leads us to conclude that our results are not an effect of extinction. 



6 



Gonzalez et al.: Detection of the inner Galactic bar traced by the VVV survey 



l = + 10° i =+5° i = -5° l = -10° 




2 1 0-1-2 
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Fig. 3. Position of the Galactic bar with respect to the Sun as measured with the RC method assum- 
ing a mean RC magnitude of = -1.55. Red and black filled circles show the results for VVV 
data at latitudes b = -1° and b = +1°, respectively. SoUd lines show the distance spread along each 
line of sight correcting for an intrinsic bulge dispersion of 0. 17 mag and photometric errors. Black 
open squares are the results from N05 at b - +1. Dashed hnes represent the o bserved orientation 



angles for true orientations of 15°, 30°, and 45° following 
the hnes of sight for longitudes I = ±5 and / ± 10. 



Stanek et al 



(11994 . Dotted hnes show 



At larger Galactic latitudes (\b\ > 3°), the flattening of the bar within |/| < 4° is not observed. 
This, combined with the symmetric change in slope, leads us to suggest that we are detecting an 
inner second bar in the Milky Way. We cannot however exclude that the observed flattening is due 
to a particular change in the density distribution of RC stars in the inner bulge. 



5. Conclusions 

We have used the VVV data at b = +1° to build the bulge luminosity functions along different 
lines of sight, and to measure the mean dereddened Kg^ magnitude, which can be used as a distance 
indicator. The RC is clearly detected in all fields, and within \l\ < 5° a second component is detected 
in the K-band luminosity function. Given that the average magnitude shift of this second peak does 
not follow the primary RC, it could be an independent structure detected at ~ 11 .2 kpc from the 
Sun e.g. a spiral arm behind the bulge. 

The main RC traces the mean orientation for the Galactic bar, which produces a magnitude 
variation from K,^ -13.4 (~ 9.6 kpc) at I = -10 to -12.4 (~ 6.2 kpc) at I = +10°. However, 
a different position angle is observed for the central ~1 kpc regions with longitudes -4° < I < 
4°. This is in excellent agreement with the results of N05, which show the same change in the 
orientation of the bar at positive latitudes b = +1°. These results provide additional evidence of a 
possible secondary inner bar with a semi-major axis of ~ 500 pc, that is symmetric to the Galactic 
plane and has an orientation angle that differs from the large scale bar. Our measurements are 
consistent with an angle of the inner bar larger than 45°, although in these regions we measure the 
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superposition of both the inner and the large scale bar and therefore detailed models are required 
to constrain the geometry and structure in the iimer bulge. 
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